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The process of dynamical reduction of the Vlasov-Maxwell equations leads to the introduction of 
classical zitterbewegung effects in reduced plasma dynamics. These effects manifest themselves in the 
form of an asymmetric canonical energy-momentum tensor involving the decoupling of the reduced 
kinetic momentum p from the reduced velocity u (i.e., u X p ^ 0) as well as reduced polarization 
and magnetization effects. The reduced intrinsic torque generated by the antisymmetric part of the 
canonical energy-momentum tensor, which is calculated from the reduced ponderomotive potential, 
acts as the source for the intrinsic (spin) angular momentum. 



In 1930, Schroedinger [l[ provided a simple interpre- 
tation for the electron's spin angular momentum based 
on Dirac's theory of the electron 01. According to 
Schroedinger's zitterbewegung model [3H5||, the intrinsic 
spin angular momentum S^" = ha^/2 of the electron 
can be interpreted as the result of the rapid trembling 
motion of the electron. The main features of the La- 
grangian formulation of zitterbewegung can be summa- 
rized by: (I) the decoupling of a particle's velocity v 
from its momentum p (i.e., v X p ^ 0); (II) the exis- 
tence of polarization and magnetization effects associ- 
ated with the decoupling between fast and slow space- 
time scales; and (III) an asymmetric canonical energy- 
momentum tensor derived from the Lagrangian density. 
In the Schroedinger zitterbewegung model, these features 
are explicitly related to the electron's spin. 

Zitterbewegung effects have recently been investigated 
in many different physical systems ranging from graphene 
d,0| and photonic materials d,Q to ultracold atoms [l(J 
and cold-ion traps [lj, [t2| ■ The purpose of the present 
Letter is to show that the reduced Vlasov-Maxwell equa- 
tions, which are obtained by dynamical reduction, exhibit 
these same zitterbewegung features derived from the re- 
duced ponderomotive Hamiltonian. 

We begin with a brief review of the Schroedinger zit- 
terbewegung model of the electron, which is governed by 
the Dirac-Maxwell equations [lij 



7"D^ = 



1 ip 



7 M = — mc ip 
d^F^ = (4tt/c)J v 



where J v = (eg, J) = ec(y v ) denotes the electron four- 
current (expressed in terms of the expectation value 
(7") = (i/jIj" \ip) of the Dirac matrix 7"), = ihd^ — 
(e/cjAfj, and = ifid^ + (e/c)A M are the kinetic- 
momentum operators acting on ip and its adjoint ip, re- 
spectively, and F^ w = d^Ay — d v A^ denotes the electro- 
magnetic field. We now show that the Lagrangian for- 
mulation of the Dirac-Maxwell equations ([I} exhibit the 
three zitterbewegung features (I)-(III). 

First, the electron is described in terms of the coordi- 
nates r, its velocity v = dr/dt, and its kinetic momentum 



p. While the electron's velocity has the eigenvalues ± c, 
its expectation value (v) (expressed in terms of the grou" 
velocity of its wavepacket) is less than c in magnitude 
The motion of the electron is therefore decomposed into 
the slow (average) motion of its center of mass and the 
rapid trembling motion of its zitterbewegung position p c . 
This rapid motion yields an expectation value for the 
magnetic dipole-moment of the electron [H, [H, [H[ 



Me 



2c ^ 



dp c 
dt 



(2) 



whose magnitude \fj, e \ = eh/(2mc) is the Bohr magneton. 

Second, we discuss the electron's polarization and mag- 
netization and introduce the Gordon decomposition of 
the electron four-current Phi 



q = ft.-V-P 

j = j c + ap/a + cVxM 



(3) 



where J" = (e/2m) \{D V ) — (D )] denotes the (spinless) 
conduction four-current and the antisymmetric magneti- 
zation tensor 



2m 



(4) 



is expressed in terms of the spin matrix <T I1U = 
(i/2) (7^7" —7^7^), whose components are the polariza- 
(1) tion AP° = cP l and the magnetization M ij = e ljk cM k . 

Third, the canonical energy-momentum tensor for the 
Dirac-Maxwell equations 



F : F 

16tt 



1 

47T 



F" a F a v - (crf^D") (5) 



is derived by Noether method from the Dirac-Maxwell 
Lagrangian density C = (cy^Dn — mc 2 ) + F^ v F Vil j\§-K . 
While the electromagnetic terms in Eq. (J5J) are symmet- 
ric, the Dirac term is not. Indeed, the antisymmetric 
part T£ v = \(T» V - T^) of Eq. © is expressed as [H 



1 A 



(6) 
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where the third-rank tensor E a ^ = - E"^ : 

(7) 

is expressed in terms of the intrinsic spin angular mo- 
mentum S aP = ha afl /2. Since (c-f 1 D^> - D l ) is 
the operator-equivalent of the cross-product v l p£ in = 
(v X Pkin)fc, the antisymmetry of the Dirac-Maxwell 
canonical energy-momentum tensor is, therefore, due to 
the decoupling of its velocity and kinetic momentum. 

Belinfante [16|, [l7| recognized that the asymmetry of 
the canonical energy-momentum tensor T^ v could be 
given a physical interpretation (see also Refs. [ID, [l9j). 
based on the fact that the transformation 

= + d a R M " (8) 

leaves the energy-momentum conservation law invariant: 

= = dX\ (9) 

where the third-rank tensor R^°^ v = — R\. a A v satisfies 
^au,^ v = 0- The condition that the new energy- 
momentum tensor (|8|) is symmetric yields the following 
expression for the antisymmetric part of the canonical 
energy-momentum tensor 

/ (10) 

and, hence, the antisymmetric part T£ v of the canoni- 
cal energy-momentum tensor acts as the source for the 
antisymmetric second-rank tensor S™''"''. 

We now show that the third-rank tensor E"''"'' is con- 
nected to the spin angular-momentum tensor as in the 
Schroedinger zitterbewegung model [see Eq. ([?]) ]■ First, 
we introduce the third-rank orbital-angular-momentum 
tensor 

and the third-rank spin-angular-momentum tensor 

= B^H - daQWIM, (12) 

where QMW = R&^x" - R^> satisfies 
dlpQ^P^ = 0. Using Eq. ©, the total angular- 
momentum tensor 

jfilrt _ L p\M + S p\M = x ^ v - X -T^ (13) 
satisfies the angular-momentum conservation law 

fyjfiw = T ^_ r » = 0j (14) 

which follows from the symmetry of the new energy- 
momentum tensor ([8]). Lastly, the equation dpL^ u ^ for 
the orbital-angular-momentum tensor (lll[) is 

fyL^H = T" u -T v » = - d S f3[ ^ ] = - dpEPW. (15) 



Hence, we see that the asymmetry of the canonical 
energy-momentum tensor acts as the source of intrinsic 
(spin) angular momentum. 

We note that the antisymmetric energy-momentum 
tensor 

?T = - \ ^ aP = \ £^ r aP (16) 

can be used to define the antisymmetric torque ten- 
sor T a p = — (d a crp — 9/3<T a ), with the spatial compo- 
nents = h e " 1 ^ TQk expressed in terms of the intrinsic 
torque tensor t ^ = (0,r): 

T = vxp + ExP + BxM. (17) 

We therefore see that all three features of the zitterbewe- 
gung model are combined in Eq. (I17|) to act as the source 
of an intrinsic (spin) angular momentum. 

After a brief introduction to the Schroedinger zitter- 
bewegung model and a general discussion of angular- 
momentum conservation within a Lagrangian perspec- 
tive, we now discuss the case of the Lagrangian formula- 
tion of reduced Vlasov-Maxwell theory. The Lagrangian 
formulation of the guiding-center and oscillation-center 
Vlasov-Maxwell equations were developed over 30 years 
ago [20i - [25| . In each case, their respective energy- 
momentum conservation laws, derived by Noether 
method, exhibited an asymmetric canonical energy- 
momentum tensor. Although it was sometimes pointed 
out that this asymmetry was only apparent [23j , based on 
the knowledge that the conservation of angular momen- 
tum required a symmetric physical energy-momentum 
tensor, Dewar [2TJ ] pointed out that the asymmetry of the 
canonical energy-momentum tensor could be expressed in 
terms of an intrinsic spin angular-momentum tensor. 

The process of dynamical reduction in single-particle 
plasma dynamics and plasma kinetic theory [26[ is asso- 
ciated with the extended near-identity canonical phase- 
space transformation % : z° = z a + e {Si,z a } + 
e 2 ({S 2 , z a } + \ {Si, {Si, z a }}) + - • • , and its inverse T~ x : 
z a = P - e {5i,2»} - e 2 ({5 2 ,^} - \ {Si, {Si,z a }}) + 
• • • , generated by the scalar fields (Si , S%, ■ ■ • ) . This 
near-identity transformation introduces the reduced- 
displacement vector 

Pe = T^x - x = - e Gf + ■ • • , (18) 

defined as the difference between the push- forward T~*x 
of the particle position x and the reduced position x. 
Hence, through the reduced displacement (TT51) . the dy- 
namical reduction yields the reduced electric-dipole mo- 
ment [m 

7T e = e Pt (19) 

and the intrinsic magnetic-dipole moment 

"-hi"**?)- « 
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which is identical in form to the zitterbewegung expression 
@ for the electron's magnetic-dipole moment. 

The extended reduced Vlasov equation for the ex- 
tended reduced Vlasov distribution T is expressed as 







d € T 
~dT 



dt dz a ' 



(21) 



where 



dU 



dt 



and 



dt 



denote the canonical Hamilton equations in extended 
phase space and the reduced Vlasov distribution is de- 
fined as T(z) = cS[w — H(x, p, t)] F(x, p, t), with the re- 
duced extended Hamiltonian satisfying the physical con- 
straint H — H — w = 0, where the reduced Hamiltonian 



H = H 



dSi 
~dt 



dS 2 
dt 



I J dSx 

2V 1 ' ~dT 



(23) 



is expressed in a form that is the classical equivalent of 
the Dirac Hamiltonian derived by the Foldy-Wouthuysen 
(FW) transformation flij . 

The dynamical reduction associated with the phase- 
space transformation 71 introduces polarization and mag- 
netization effects into the Maxwell equations, which are 
transformed into the macroscopic (reduced) Maxwell 
equations [25l [2(| 



1 <9D 47r - 

V-D = 47rg and VxH — = — J, (24) 

c dt c 

where the microscopic electric and magnetic fields E and 
B are replaced by the macroscopic fields D = E + 4ir P e 
and H = B — 47rM e , where P e and M e are the reduced 
polarization and magnetization. We note that the dy- 
namical reduction associated with the phase-space trans- 
formation 71 has introduced the following expressions for 
the charge and current densities: 



Q 
J 



Q- V-P e 
J + dPJdt 



:VxM f 



(25) 



which are of course similar to the Gordon decomposition 
§3§ observed in the Dirac model. 

Lastly, by using the reduced electric-dipole and 
magnetic-dipole moments (fT9)) - (f2"0)) , we construct explicit 
expressions for the reduced polarization 

Pe = E / ^Fd 3 p, (26) 
and the reduced magnetization 

+ MR Pi) 

which combines the intrinsic magnetic-dipole contribu- 
tion (|20[) and the moving electric-dipole contribution. 



We now show that the reduced Vlasov-Maxwell equa- 
tions (|2"Tj) and (|24p can be derived from the reduced vari- 
ational principle J d^x 5C — 0, where the reduced La- 
grangian density is 



£( X ) = Yl 
1 



d A p T[x, p) H(x, p; A, F) 



+ —F(x):F(x). 
Ion 



(28) 



Note that, as a result of the dynamical reduction of the 
Vlasov equation, the reduced Hamiltonian (|2T)1) is ex- 
pressed in terms of canonical energy-momentum coor- 
dinates as 



H(%p,t;A, F) = -L|p-(e/c)A| i 
2m 

- A; F 

c 



4> 



(29) 



where the reduced ponderomotive potential ty e depends 
explicitly on the field tensor F M „ , which once again is the 
classical equivalent of the Dirac Hamiltonian after the 
FW transformation. From these field dependences, we 
define the reduced four-current density 



J" = (eg; J) 

where d e x°/dt = c and 



d p J 7 



d e x» 
dt ' 



d £ x dU 1 _ e A , <9* e 
-jr = -q^t = -(p--A) + —. 
dt op m c op 



(30) 



(31) 



The reduced antisymmetric polarization-magnetization 
tensor 24], on the other hand, yields the reduced polar- 
ization = PI and the reduced magnetization K l J = 
e ljk M e k , where 



(Pe, Me) 



E 



dpF 1 flE» 9B 



(32) 



Hence, polarization and magnetization are explicitly 
evaluated in terms of the ponderomotive Hamiltonian ^> c . 

By applying the Nother method, we obtain the reduced 
canonical energy-momentum tensor 



JflV = _ 
4-7T 



9 



F : F - F 



47rKr F 



E/^(p 



-A") r, 



(33) 



which naturally includes the polarization and magneti- 
zation (1321) and is manifestly asymmetric. The antisym- 
metric part of the canonical energy-momentum tensor 
T 1 ^ = | e^ k Tk is expressed in terms of the reduced (in- 
trinsic) ponderomotive torque 

~8$ r 



E 



<9E 



X E 



~dB 



x B 



<9* e e 



Fd 3 p. (34) 



4 



The three zitterbewegung features, which appear explic- 
itly in the expression for the reduced ponderomotivc 
torque, are derived from the reduced ponderomotive po- 
tential 9 e associated with the process of dynamical re- 
duction. 

As an application of the ponderomotive torque (|34p . 
we consider the gyrocenter ponderomotive Hamiltonian 

M 



9 e 



mc 
2Bl 



E, 



P||bo 



x B, 



IB 



2Bn 



-, (35) 



derived in the zero-Larmor-radius (ZLR) limit from 
the second-order gyrocenter Hamiltonian [281 ] ■ where a 
strongly magnetized background plasma is perturbed by 
the low-frequency electromagnetic fields E_l = -Vi$ 
and B^ = Vj_-A|i X bo, which are perpendicular to the 
background magnetic field Bo = So bo- From the gyro- 
center ponderomotive Hamiltonian (|35p . we obtain the 
reduced gyrocenter polarization and magnetization 



d9 e 
9E ± 

d9 e 
where /i e 



mc 



E 



P||bo 



mc 



x B i 



7T e X 



Piibo 



mc 



7T e , (36) 



(37) 



/iB^/Bo denotes the reduced in- 



trinsic gyrocenter magnetization. These definitions 
lead to the identity d9 e /d~E± X E ± = Ej_ X iv e = 
— d9 e /d'B± X B^, which means that the electromag- 
netic contribution to the gyrokinetic intrinsic torque van- 
ishes in the ZLR limit. The kinetic part of the gyrokinetic 
intrinsic torque, on the other hand, involves the pondero- 
motivc velocity 



d9 e 
dp 



^_ [ E £bo Pj| B_l 

Bo 1 B m B Q 



99 f 



dp 



(38) 



which vanishes only in the electrostatic limit (Bj_ = 0). 
Lastly, the intrinsic gyrokinetic torque 



gy 



E 



89, 



dp 



p X b F d 3 p 



(39) 



involves the magnitude of the ponderomotive velocity 
and the perpendicular component of the gyrocen- 



ter kinetic momentum [p — (e/c) A\\ bo] X bo = p X bo- 

In summary, we have shown how classical zitterbe- 
wegung effects manifest themselves in reduced Vlasov- 
Maxwell theory through the reduced ponderomotivc 
Hamiltonian. Here, the ponderomotive torque (f34|) acts 
as the source for the intrinsic spin angular momentum. 
Ongoing work includes the investigation of classical zit- 
terbewegung effects in gyrofluid and gyrokinetic models. 
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